Collective modes of the d-density wave state and its relevance to high-T c cuprates 
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We calculate the collective mode spectra of the d-density wave state, proposed to describe the anomalous 
phenomenology of the cuprates in the pseudogap regime. Even though the state breaks translational symmetry 
by a lattice spacing and is described by a particle-hole singlet order parameter condensing at the wave vector 
q = Q = (tt, 7r), interestingly, we find that the amplitude collective mode spectrum can have inelastic peaks at 
both q = (0, 0) and q — Q— (7r, 7r). In general, the spectra is non-universal, and, depending on the microscopic 
parameters, can have one or two peaks in the two-dimensional Brillouin zone, signifying confluence of two kinds 
of magnetic excitations. Our theory should be important in interpreting inelastic neutron scattering experiments 
in high-T c cuprates, where such multiple magnetic excitations are implicated. 
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Ever since the discovery of pseudogap in high temperature 
superconductors it has been a profound mystery Q]. To this 
day, its origin is vigorously debated. In one view, pseudogap 
is a remnant of the d-wave superconducting gap that defines 
a crossover temperature T* in the phase diagram. The other 
view argues for a broken symmetry at T*. The precise na- 
ture of the broken symmetry is debated, however (321 . Here, 
we shall assume that much of the phenomenology associated 
with the pseudogap can be described in a unified manner by 
the single assumption of a spin singlet d x 2_ y 2 density wave 
(DDW) [3 ]. The purpose of the present Letter is to deduce 
certain observable consequences of this assumption, namely 
multiple collective modes resulting from this ordered state, 
which can possibly arise in inelastic neutron scattering exper- 
iments, as in Refs. 0[6). As we shall see, these turn out to 
be unusual and were missed in a very early paper [ 7 ] on this 
subject. As to elastic signature of singlet DDW, two neutron 
scattering measurements seem to provide some evidence for 
it 1 8 ]. However, these measurements have not been confirmed 
by further independent experiments. 

The singlet DDW, as originally envisioned, consists of cir- 
culating currents (9|, alternating between clockwise and anti- 
clockwise directions in the neighboring plaquettes of an un- 
derlying square lattice in two-dimensional (2D) CuC>2 planes. 
Viewed from this perspective, the state can be constructed as 
a juxtaposition of two kinds of current conserving vertices, as 
shown in vertices 5 and 6 of Fig. [T] The ordered state breaks 
time reversal, rotation by tt/2, parity, and translational sym- 
metry by one lattice spacing, but the product of any two of 
these symmetries is preserved. Superficially, the statistical 
mechanics belongs to the Ising universality class 0, and in 
mean field theory the order disappears when the magnitude 
of the bond currents vanishes with increase of temperature. 
However, fluctuations can reverse an arrow if it is possible to 
do so in a current conserving manner. This can be done by 
enlarging the configuration space by adding four additional 
vertices shown as vertices 1, 2, 3, and 4 in Fig.[T][ 10 ]. The ver- 
tices corresponding to sources and sinks should have a large 
negative chemical potential. The model then belongs to the 
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FIG. 1. The six possible current vertices in the 6-vertex model. 
The vertices (5) and (6) are the AF vertices which lead to the DDW 
phase with local orbital moments, while the rest of the vertices lead 
to longer range current fluctuations. 



universality class of the classical six-vertex model [10]. The 
quantum six- vertex model and its quantum criticality has also 
been previously discussed ifTTTl . 

The particle-hole spin- singlet DDW order parameter is 



( C k,c C k+Q,/3) = iA Q fk5 i 



a/3, 



(1) 



where c, are electron annihilation and creation operators, 
and a, f3 are the spin indices; Aq is the magnitude of the or- 
der parameter and the form factor f^ = (cosk x — cosk y ); 
we set the lattice constant to unity throughout. Because the 
order parameter condenses at the wave vector Q — (7r,7r), 
one would have naively expected that the amplitude collective 
mode, the analog of the U(l) Anderson-Higgs mode, but in 
the particle-hole channel, would be peaked at the same wave 
vector. 

An important result of our paper is that the spectra of the 
amplitude fluctuations of the DDW state are not confined to 
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only the ordering wave vector Q, but can have finite-frequency 
peaks at both q = (0, 0) = and q = Q, as well as con- 
siderable spectral weight over a substantial region of the mo- 
mentum space. The emergence of a q = peak, even with 
a mean field state that breaks the lattice translation symme- 
try such as the twofold commensurate DDW, indicates q = 
fluctuations hidden within the DDW state. We find that these 
fluctuations reflect an orbital ferromagnetic (OF) phase in the 
six-vertex generalization of the model, as mentioned above. 
In its phase diagram one finds OF, as well as orbital antiferro- 
magnetic (DDW) phases ifTOl . The existence of such striking 
q = fluctuations, even in the absence of an OF order param- 
eter, can be important in interpreting inelastic neutron scat- 
tering experiments [5]. However, because the response is at 
higher frequencies, they are subject to considerable degree of 
non-universality and serves as a warning that if our theory is 
correct, experimental signatures should not be unique across 
materials, as well as within a given material, as a function of 
doping, frequency, etc. 

To calculate the DDW collective mode spectrum we start 
with the electronic Hamiltonian which has the form, 

H = - fi)cl a c k , a + g- 1 A g , (2) 

k,a q 

where 

4z = fZ) fk+{q-Q)/2C ] k C k+q , (3) 
k 

Here we adopt a commonly used dispersion e k based on local 
density approximation [12], which is 

e k = — 2t(cos k x + cos k y ) + 4t' cos k x cos k y 
— 2t" (cos 2k x + cos 2k y ). 

For subsequent notational simplicity, we define the functions, 
Sk = \{?k + e fc+ Q), and e k = \{e k - e k+Q ). The band pa- 
rameters are chosen to be t = 0.15 eV, t' = 0.3£, and 
t" — 0.5t'. The difference with the conventional LDA band 
structure is a rough renormalization of t (from 0.38 eV to 0.15 
eV), which is supported by experiments involving angle re- 
solved photoemission spectroscopy fT3l . Note that in Eq. [3] 
as elsewhere, we shall drop the spin-index a, as it will play no 
important role. 

For an analysis of the collective modes of the DDW state at 
a general q, we have generalized the gap parameter in Eq. [3] 
In this form, A J = A _ q , for an arbitrary q, which implies 
that the gap parameter is purely real in the real space. This 
constraint is necessary to ensure that no gapless phase mode 
is generated in the collective mode spectrum of the twofold 
commensurate DDW state which breaks only discrete sym- 
metries Q. From Eq.jijit follows that at q = Q, Aq repre- 
sents the conventional DDW gap parameter shown in Fig. [T] 
in its vertex representation, while q — leads to a gap pa- 
rameter A which represents uniform current flow along the 
+x and +y directions. Choosing q = Q = Q + Q, where 
Q = (7r, — 7r), we find currents flowing along the +x and — y 



directions. These, together with reversal of currents, give the 
remaining vertices in Fig.[l] Finally, at q = Q, A q allows a 
gap parameter that breaks local current conservation with sites 
having a current source at one vertex and a current sink in the 
neighboring vertex. However, in this paper we will choose 
the fugacity of these current-conservation violating vertices, 
which are controlled by the coupling constant g in front of the 
form factors, to vanish [10]. Thus the operator A q represents 
the full set of current vertices of the six- vertex model. We find 
that only the DDW gap parameter Aq develops a mean-field 
expectation value in the saddle-point solution, though we will 
find fluctuations from all the other vertices as well. 

Anticipating a twofold commensurate DDW order with 
the order parameter given by Eq. ([I]), we first fold the full 
2D Brillouin zone (BZ) to the reduced BZ (RBZ). The re- 
duced zone is defined in terms of the rotated coordinates, 
{k x + k y )/\^2 = k' x , (k x — k y )/y/2 = k' yJ so that in the 
RBZ k' x , k f y e [-?r/V% k/V2}. Note that the Dirac points in 
the spectrum occur at (k' x , k' ) = (n/y/2, 0) and (k' x , k' ) = 
(0, ir/y/2). The vector Q is at (V^tt, 0), while Q = (tt, -tt) 
in the original basis is now at (0, ^tt). 

To facilitate our discussion, we introduce the spinor nota- 
tion & k = (<4,<4 + q). Defining a BZ periodic function u k 

which is 1 inside the RBZ and zero outside, A q may be writ- 
ten as (ai, i = 1, 2, 3 are the conventional Pauli matrices) 

= ^2 U kfk+(q-Q)/2^{[Uk+qCr3 + iu k+q+Q a 2 ]^ k+q 
k 

The above expression for A q is defined for q in the full BZ. In 
the RBZ this corresponds to four different gap parameters. It 
is convenient to split A q as 

A g = ^ U k U k j rq f k _Q/2+ q /2^\v^k+q 
k 

\+Q = ^2 U k U k+qfk+q/2^\v2^k+q 
k 

Ag+Q = ^U k U k+q f k _Q /2+q/2 ^\(J^ k+q 
k 

Ag+Q + Q = ^U k U k +qf k+q / 2 +Q/2 1 ®{cr2 1 ®k+q- (6) 

k 

where Q = (Q + Q) and we have dropped the imaginary i 
in A q +Q and A q+ Q + Q. After this transformation, the vec- 
tor q only takes values which are the differences of the wave- 
vectors in the RBZ. Physically, these 4 operators (with posi- 
tive and negative signs) now correspond to the 8 current order- 
ings of the eight-vertex model. As discussed above, we will 
choose the g associated with the non-current conserving ver- 
tices to vanish, thus leaving us with a six-vertex model. The 
phase transition of the six-vertex model does not have a sin- 
gular specific heat at its transition, whereas the eight vertex 
model does; this is in accord with experiments in high tem- 
perature superconductors. 

For compactness of notation we will write A^ = 
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A g+ 5+ a Q> where a = 0, 1 and b = 0, Q. With this nota- 
tion, the hermiticity condition on A q translates into the condi- 
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The collective modes are calculated by considering the 
terms second order in A q ^ in Eq.[9] which have the form 



. Furthermore it is convenient to write 
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£ fc *ilAL/2*^ where #2+*/2 are 2 x 2 ^ (2) ~ ^ ? 
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matrix- valued structure factors. The Hamiltonian in the spinor 
notation is 

H = ^*l((-M+e fe )^+ e ^ 3 )* fe +5- 1 E [ A M tA i:<! 

(7) 

where in the sums over a and b, a = 1 and 6 = Q are 
dropped because these vertices do not conserve currents. Con- 
sider now the imaginary-time effective action and decouple 
the four-fermion term S g a & a£^A^ in Eq. 7 using the 
Hubbard- Stratonovich transformation. This leads to the ef- 
fective action, 

S = ^ ^k,co((~ iuj ~ V + £ k) cr + £kCT3)^k,oj 



k.oj 



2 * _1 E A £t A(a) 



b,q,uj 



-1 A *0) a (a) 

q,oj,a,b 



(8) 



where the Hubbard- Stratonovich fields satisfy the constraint 
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± b _ q _ UJ = Aj*^. Expressing the A operators in Eq. 
in terms of the fermion spinors \I/'s and then performing tii 
Grassmannian path integral over the quadratic terms leads to 
the effective action, 

o _ -l A *(o) a (a) 

q,oj,a,b 

-TiKM 0>fcl , Wl i fcl>fc2 i Wl ^+<^)], (9) 

where 

M 0iki Lj = (-iu - V + £k)<ro + Ck(?3 - y^Pb a l A b a pi ( 10 ) 
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The mean-field equation for the DDW state is obtained by 
setting the terms which are linear in SM in Eq.|9]to zero. This 
leads to the equation, 



1 fk r /Zk + Ek-ii e k -E k 



2T 



2T 



(12) 



for the order parameter Aq; E k = y + fk^Q- and n^i^) 
is the Fermi-function. The chemical potential fi is determined 
from the hole-doping x. There are no static saddle point so- 
lutions corresponding to (A q ) and (A^+q) for the range of 
parameters considered. However, interestingly, the collective 
mode spectrum for the DDW state will contain a component 
also near q ~ 0, which can be interpreted as orbital ferromag- 
netic fluctuations of the bond currents. 



TrK^ 1+ J^pg i+ ^ 

a, b 

a, 6 

a;, a, 6 



(13) 



Performing the Matsubara summations and analytically con- 
tinuing to real frequency leads to 

a,b,a f ,b' 

where /7^' a ^ (g, a;) are given by 

Uf£, ,a \q,u) = ^2 A m , n (/ci,^i + g;o;) 

/ci ,m,n 



(15) 



and 



A m , n (fci,fci + = 

n F (mE kl + g fcl - /i) - n F (nE kl + q + g fcl+g - /i) 
2 (cj - m£?fei + nE kl + q - e kl + e^+g + iS) 

(16) 

with A m>fe = ^(E k ^me k a 3 -mf k AQa 2 ). Fortherestof 
the paper we will focus on the spectra at the special q points 
q = Q and q = where the U coefficients take the form 



<q\^)= E 



k:E k >\fi-£ k \ 



E k (w*-4E%) 



(17) 



and 

T (a,a) 



k:E k >\fi-£ k \ 



E k {^-AEl) 



(18) 



respectively at T = 0. In this limit, we find that all terms, 
where b ^ b' or a ^ a', vanish because of the different sym- 
metries of the form factors under k x — » —k x and k y —k y . 

The measurement of dissipation from any probe that cou- 
ples to the fields A q and A g+ Q would be of the form 



(a,a) , 
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with 



V(q,oj) 



(20) 



where (a, b) and (a', b') are treated as matrix indices. Thus the 
peaks in the imaginary part of the spectrum near q = and Q 



are given by the zeroes of the functions in Eq.[T7]and Eq.[18 
The a = (corresponding to the conventional DDW order 



parameter) propagator near Q, given by Eq. 17 can be easily 
seen to vanish when uj ~ 2Aq because of the numerator. On 
the other hand, low uj zeroes in the other propagators only 
emerge when \i is such that the denominator (uj 2 — 4£|) can 
vanish while satisfying the condition Ek > |/x — Sk\- This 
occurs in the presence of a non-zero next nearest neighbor 
coupling t f , which breaks particle-hole symmetry in the band 
structure. 

There is no mean field OF order coexisting with the twofold 
commensurate DDW order in the regime of doping studied 
here. In contrast, from Fig. [2] it is clear that in the under- 
doped regime there is a noticeable finite-frequency peak at 
q = coexisting with a finite frequency peak at Q. The in- 
tensity at q = goes down with increasing x, but the ampli- 
tudes strongly depend on the microscopic parameters, which 
implies that the amplitude fluctuation spectra can have finite 
frequency peaks at both wave vectors q = and q = Q, or 
only a single peak. Such non-universality indicates that dif- 
ferent families of cuprates, or even samples at different hole 
doping within the same family of cuprates, may have differ- 
ent fluctuation spectra. The non-universality of the collective 
modes and the coexistence of finite frequency peaks at mul- 
tiple wave vectors are the central aspects of our work. These 
results are direct consequences of the vertex model framework 
of the fluctuations for the DDW state and may have important 
consequences for the inelastic neutron scattering experiments, 
which are left for future study. 

In summary, we have calculated the collective mode spec- 
tra for the commensurate DDW state, which was proposed 
earlier to explain the phenomenology of the pseudogap phase 
of the high-T c cuprates. To go beyond the fluctuations at Q, 
the DDW gap parameter was generalized as in Eq.[3] The per- 
sistence of a hidden, q = OF fluctuations even when the 
ground state is the twofold commensurate DDW state is nat- 
ural from the vertex model description of the DDW state, In 
principle, one can also expect fluctuations associated with the 
q = nematic state, which, as shown in Ref. fl4l . is compara- 
ble in energy to the DDW state. This nematic order parameter 
Anem = Ylk fk c \ c k is a close cousin of the order parame- 
ter discussed here. Finally, because the orbital magnetic mo- 
ment contribution to inelastic neutron scattering is expected 
to be enhanced at long wavelengths, the relatively small peak 
at q = can produce scattering peaks comparable to that at 
Q. However, the detailed kinematics addressing quantitatively 
the neutron scattering signatures of the collective modes are 
beyond the scope of the present paper. In any event, because 
of the non-universality of the inelastic spectra, as mentioned 
above, such an effort will not be particularly meaningful at 
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FIG. 2. (Color online) Collective mode response of the commen- 
surate DDW state as a function of uo/t for q = Q (solid curve) and 
q = (dashed curve) for the band parameters given in the text. In the 
panels (a)-(c), A Q = 0.3t, and x = 0.06 (a), x = 0.1 (b), x = 0.14 
(c), x — 0.18 (d). The emergence of an inelastic q = fluctua- 
tion peak with underdoping is a signature of hidden underlying OF 
fluctuations co-existing with the commensurate DDW fluctuations 
centered at q = Q. The vertical axes correspond to arbitrary units, 
but are the same for all panels.The panels (e) and (f) correspond to 
x — 0.06 and Aq = 0.4t and Aq = 0.5t respectively. 
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